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Abstract 


A parametric study of the effects of test-fixture-induced initial prestress 
and elastic edge restraints on the prebuckling and buckling responses of a 
compression-loaded, quasi-isotropic curved panel is presented. The numerical 
results were obtained by using a geometrically nonlinear finite element analysis 
code with high-fidelity models. The results presented show that a wide range of 
prebuckling and buckling behavior can be obtained by varying parameters that 
represent circumferential loaded-edge restraint and rotational unloaded-edge 
restraint provided by a test fixture and that represent the mismatch in specimen 
and test-fixture radii of curvature. For a certain range of parameters, the panels 
exhibit substantial nonlinear prebuckling deformations that yield buckling loads 
nearly twice the corresponding buckling load predicted by a traditional linear 
bifurcation buckling analysis for shallow curved panels. In contrast, the results 
show another range of parameters exist for which the nonlinear prebuckling 
deformations either do not exist or are relatively benign, and the panels exhibit 
buckling loads that are nearly equal to the corresponding linear bifurcation 
buckling load. Overall, the results should be of particular interest to scientists, 
engineers, and designers involved in simulating flight-hardware boundary 
conditions in structural verification and certification tests, involved in validating 
structural analysis tools, and interested in tailoring buckling performance. 


Introduction 

Curved panels are a common structural element found in many aerospace vehicles. With the need 
to reduce structural weight of aircraft and to exploit stiffness tailoring of fiber-reinforced, laminated- 
composite materials to enhance performance, there has come the need to understand the behavioral 
characteristics of thin-walled curved panels that may possess high degrees of material orthotropy and 
anisotropy. Orthotropy and anisotropy are manifested physically by couplings between deformation 
modes such as extension, bending, and twisting. These coupled deformations may enable the design of 
special-purpose, high-performance components that meet strict structural requirements. Because these 
structures are typically thin walled and are often subjected to compression loads, predicting their 
prebuckling and buckling behaviors is of fundamental importance. A necessary step in developing this 
predictive capability is the development of high-fidelity analysis capabilities and the corresponding 
validation procedures. 

Significant advances in high-fidelity nonlinear finite-element analysis methods and validation- 
experiment requirements for buckling of thin-walled, laminated-composite cylindrical shells have been 
made in recent years (e.g., see ref. 1). These analysis methods are mature and lend themselves to 
compression-loaded curved panels, but the validation of these methods for curved panels is much less 
mature. The primary reasons for this lack of maturity, compared to that for compression-loaded 
cylindrical shells, are that laboratory testing of curved panels requires more complex fixtures to facilitate 
load introduction and edge support. Furthermore, the effects of the test fixture on the load introduction 
and edge support conditions can be difficult to characterize. These issues are significant because careful 
attention to the implementation and characterization of load introduction and edge support conditions 
used for validation tests is necessary to achieve relevancy and usefulness of a validation experiment in the 
development of practical structural design data. In particular, the test results may not be relevant if the 
test configuration does not accurately represent the structural configuration in question. Similarly, the test 



data may be of limited use for model validation purposes if the load introduction and edge support 
conditions are not accurately characterized and thus are not properly accounted for in the analysis. 

Many theoretical and experimental studies that are relevant to validation of high-fidelity analysis 
methods for bucking of compression-loaded curved panels have been conducted during the last thirty 
years (refs. 2-8). In particular, these studies have examined the effects of panel curvature, initial 
geometric imperfections, boundary conditions, and laminate stacking sequence on the buckling behavior. 
Many of these studies have shown significant discrepancies between the experimentally determined 
buckling loads and the corresponding analytical buckling predictions (refs. 5, 6 and 8). These 
discrepancies have been attributed to various effects such as improper simulation of boundary conditions, 
initial geometric imperfections, and initial prestress caused by mismatches between panel and test-fixture 
geometries. For example, experimental results obtained by Hui (ref.4) for buckling of symmetric cross- 
ply laminates show that compression-loaded curved panels can exhibit buckling loads greater than the 
corresponding classical, linear-bifurcation buckling load if an initial geometric imperfection in the form 
of an outward bulge is present in the panel. This result is in complete agreement with the unstable, 
asymmetric bifurcation behavior predicted by even simple analytical methods for compression-loaded 
curved panels (ref. 9). 

Khot and Bauld (ref. 5) conducted an experimental study of the buckling response of 
compression-loaded curved panels made from specially orthotropic laminates that is particularly relevant 
to the development of validation experiments. In this study, the experimentally measured buckling loads 
were found to be bounded by the predicted buckling loads for the case where the circumferential 
displacements on the curved loaded edges are unrestrained (free circumferential expansion) and for the 
case where these displacements are fully restrained. The bounding of the experimental buckling loads by 
the analytical predictions for these two extreme, idealized boundary conditions was attributed to the 
presence of circumferential slippage of the panel edges in the test fixture, which was observed during the 
tests. It was also noted in this study that some of the geometrically imperfect panel specimens had to be 
forced slightly into the more geometrically precise test fixture, which produced initial stress and 
deformation states prior to compression loading. It was speculated that perhaps these test-fixture-induced 
states could have an influence on the buckling resistance that is comparable to that of an initial geometric 
imperfection. 

Wilkens (ref. 6) conducted extensive buckling tests of compression-loaded, laminated-composite 
curved panels and suggested that friction between the test fixture and the specimen can affect the panel 
response. Wilkins minimized the friction that develops between the test fixture and the specimen by 
bordering the specimen edges with 0.003-in-thick Teflon tape. In addition, Wilkins showed that test- 
fixture edge supports that are intended to simulate simply supported boundary conditions on the unloaded 
edges actually provide a slight rotational constraint. For the curved panels he considered, the maximum 
value of the restraining moment was found to be approximately 1 0 in- lb/in. along the length of the panel 
edge. However, this restraining moment was determined to have a negligible effect on the magnitude of 
the buckling load of the panels considered. 

A recent numerical and experimental study by Hilburger et. al (ref. 8) also investigated the effects 
of friction-induced circumferential edge restraint on the prebuckling and buckling response of 
compression-loaded, quasi-isotropic curved panels. Specifically, experimental results were obtained that 
revealed a significant amount of nonlinear prebuckling deformation. In addition, the experimental 
buckling loads were nearly twice the classical linear-bifurcation buckling load predicted by using a finite- 
element analysis. The experimental buckling load was found not to be bounded by the corresponding 
predicted buckling loads for the case where the circumferential displacements on the curved loaded edges 
are unrestrained (free circumferential expansion) and for the case where these displacements are fully 
restrained. This inability to bound the results for these two extreme boundary conditions led to the 
development and use of a finite-element model that simulates test-fixture-induced frictional forces by 
using circumferentially aligned elastic springs distributed along the curved, loaded edges of the panel. 

2 



The numerical results obtained with this enhanced finite-element model indicate that a panel with 
circumferential edge restraint provided by the elastic springs can exhibit significant nonlinear prebuckling 
deformations that are caused by a nonlinear coupling between the circumferential and out-of-plane 
displacements near the loaded edges. Additionally, the nonlinear prebuckling deformations result in a net 
radially outward bulge of the panel prior to buckling that increases the buckling resistance in the same 
manner as reported by Hui (ref. 4) for initial geometric imperfections. Moreover, the results obtained by 
using the springs at the loaded edges appear to more accurately represent the experimentally measured 
results and suggest that some amount of circumferential restraint of the loaded edges was present in the 
experimental setup. 

The results presented in the literature reviewed herein indicate several factors related to load 
introduction, boundary conditions, and geomtery that must be well understood in order to develop 
experiments that are adequate for validating analysis methods for laminated-composite curved panels. 
Three factors considered in this paper are the effects of circumferential edge restraint, the effects of 
rotational edge restraint, and the effects of a geometric mismatch between the panel and test- fixture at the 
support locations. Thus, the objective of the present study is to quantify these three effects, and their 
interactions, for the baseline case of a compression-loaded, quasi-isotropic composite curved panel by 
using high-fidelity nonlinear finite-element analyses. Toward this goal, results of a parametric study are 
presented that encompass a wide range of circumferential and rotational edge-restraint conditions and 
initial stress and deformation states associated with forcing a geometrically imperfect panel specimen into 
a more geometrically precise test fixture, prior to application of the destabilizing compression load. First, 
the finite-element model and analysis method are described. Then, results that illustrate the individual 
contributions of each of the three effects on the prebuckling and buckling response are presented and 
discussed. Lastly, results that show the nonlinear interaction of the three effects are presented. 


Finite-Element Model and Analysis Method 

The panels considered in this study were analyzed with the STAGS (STructural Analysis of 
General Shells) nonlinear shell analysis code (ref. 10). STAGS is a finite-element code designed for the 
static and dynamic analysis of general shells, and includes the capability to model the effects of geometric 
nonlinearities. A typical finite-element model of a curved shell is shown in figure 1. Points making up 
the panel mid-surface are located by an x-y-z curvilinear coordinate system, with the origin at the center 
of the panel. The mid-surface axial, circumferential, and out-of-plane displacements are defined as u(x, 
y), v(x, y), and w(x, y), respectively, as shown in figure 1. Here, the term "out-of-plane" refers to the 
displacement that is normal to the tangent plane at a given point of panel mid-surface. For the case 
studied, the nominal panel radius R is 60 in. and the panel length L and arc width W are 14.75 in. and 14.5 
in., respectively. The panel is modeled as a 24-ply-thick [±45/0/90] 3s quasi-isotropic graphite-epoxy 
laminate with a nominal ply thickness of 0.005 in. and a total, nominal panel thickness of 0.12 in. The 
lamina material properties are as follows: longitudinal modulus E\ = 18.5 Msi, transverse modulus E 2 = 
1.64 Msi, inplane shear modulus Gu = 0.87 Msi, and major Poisson’s ratio v [2 = 0.304. 

Idealizations of the test-fixture support conditions were used in the finite-element model, referred 
to herein as the nominal boundary conditions, and are defined as follows. On the clamped loaded edges, 
the out-of-plane displacement w, is set equal to zero in the boundary regions of the finite-element model 
that extend 3/8 in. from both loaded ends of a panel (see figure 1). The compression load is introduced 
into the panel by applying a uniform end shortening displacement A to one end of the model while 
holding the other end of the panel fixed on the boundary; that is, u(-L/2, y) = A and u(L/2, y) = 0. The 
simply supported boundaries on the unloaded edges (knife-edge supports) were simulated by setting the 
out-of-plane displacement, w, equal to zero on a line 3/16 in. from each unloaded edge of a panel. The 
finite-element model uses four-noded, flat facet-type elements that are based on the kinematic hypothesis 
of Kirchoff-Love thin-shell theory and the nonlinear Lagrangian strain tensor. The element nodes include 
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three translational degrees of freedom ( u , v, w) and three rotational degrees of freedom. The rotational 
degrees of freedom are denoted by r x , r y , and r- and correspond to rotations about the x, y, and z axes, 
respectively. Large rotations are accounted for by using a co-rotational algorithm. 

The STAGS finite-element code uses both the modified and full Newton methods for its 
nonlinear solution algorithms. The Riks pseudo arc-length path- following method (ref. 10) is used to 
continue a solution past the limit points of a nonlinear response. In this computational approach, the 
incrementally applied loading parameter is replaced by an arc length along the solution path, which is 
then used as the independent loading parameter. The arc-length increments are automatically adjusted by 
the program as a function of the nonlinearity exhibited in the solution. 

The analysis considers idealized, geometrically perfect panels with nominal boundary conditions 
and panels having elastic boundary restraints with initial stress and deformation states caused by a 
mismatch between the panel and test-fixture geometry. In particular, two distinct elastic boundary 
restraints are considered; that is, elastic circumferential and elastic rotational edge restraints. An elastic 
circumferential edge restraint was applied to the curved, loaded edge as a “first-order” approximation of 
the frictional forces generated between the test fixture and the panel specimen. The elastic 
circumferential edge restraint is simulated in the finite element model by adding a linear spring-stiffness 
contribution K v directly to the stiffness of the elements in the curved, loaded boundary that are indicated 
in figure 1. The magnitude of K v varied from a value of 0.0 lb/in., corresponding to unrestrained 
circumferential motion, to a relatively high value of 1.0 X 10 9 lb/in. In addition, complete circumferential 
restraint at the loaded edges, given by v = 0, was also evaluated. As a point of reference, the 
circumferential membrane stiffness (A 22 ) of the panel is 9.22 X 10 5 lb/in. 

To better understand the effects of the restraint generated by the knife-edge supports on the 
unloaded edges of the panels, an elastic rotational edge restraint was applied to each unloaded edge that is 
representative of a “first-order” approximation to the true support condition. The elastic rotational edge 
restraint is simulated in the finite-element model by a beam on both sides of the panel that has a 
prescribed torsional stiffness and values for the bending and extensional stiffness that are equal to zero. 
The magnitude of the torsional stiffness K e varied from a value of 0.0 lb-in, corresponding to free 
rotation, to a relatively high value of 1.0 X 10 9 in-lb. In addition, complete rotational restraint at the 
unloaded edges (i.e., clamped), given by r x = 0, was also considered. As a point of reference, the 
corresponding circumferential bending stiffness (D 22 ) of the panel is 1.2 X 1 0 ’ in- lb. 

Initial stress and deformation states caused by a geometric mismatch between the panel and the 
test fixture at the edge supports were also modeled. For conciseness, these states are referred to herein 
collectively as a prestress state or a prestress condition. For example, the test fixture used in ref. 8, and 
modeled in the present study, is made for a 60-in-radius curved panel, however, because of residual 
stresses from the curing process, the curved laminated-composite test specimens had a measured average 
constant radius of about 53 to 54 in. Thus, when a panel was inserted into the 60-in-radius fixture, the 
panel was forced to conform to the shape of the fixture, which generates u nk nown stresses and 
deformations in the panel. This misfit condition is modeled by first simulating a stress-free panel with an 
initial radius of R, that represents the as-built average panel radius. Then, displacements are applied to the 
edges of the panel that simulate insertion of the panel into the 60-in-radius test fixture. The deformation 
associated with inserting the panel in the fixture is accomplished by applying an out-of-plane 
displacement to all the nodes on the edges that are based on the difference between the initial, as-built 
panel radius and the test-fixture radius. After the initial stress and deformation states are applied to the 
panel, the destabilizing compression load is applied by prescribing a uniform end shortening that is 
increased monotonically from a value of zero. 

Model convergence studies and model verification studies were conducted at various stages in the 
development of the finite-element models to verify the adequacy of each model. Specifically, 
convergence studies on the mesh refinement were conducted for each panel configuration; e.g., panel 
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configurations with elastic edge restraint only, initial prestress only, and combinations of elastic edge 
restraint and initial prestress. The convergence criteria used in the study required that the buckling load 
could not change by more than one percent with additional mesh refinement. In addition, the 
convergence criteria required that the overall predicted global response characteristics could not change 
with additional mesh refinement, e.g., the load-shortening response curves and the prebuckling and 
buckling displacements. The rationale for using a convergence criteria with requirements on the buckling 
load as well as other global response characteristics is to assure convergence on the correct overall panel 
response and not just the value for the buckling load. Model verification studies were also conducted to 
insure that the model inputs were correct. For example, the displacements applied to the clamped 
boundary region of the panel that were used to simulate the mismatch between the panel and the fixture 
geometry were verified by comparing the displacements reported in the finite-element model output file 
with the displacements predicted by basic analytical geometry. 


Results and Discussion 

Numerical results are presented in this section for 60-in-radius, quasi-isotropic panels that are 
subjected to uniaxial compression loads. First, results are presented that show the effects of elastic 
circumferential edge restraints applied to the curved, loaded edges of the panels. Similar results that have 
been presented in ref. 8 are included herein for completeness. Then, results are presented that show the 
effects of elastic rotational restraints applied to the straight, unloaded edges of the panels. Next, results 
are presented that show the effects of the initial prestress condition that is caused by the panel-test-fixture 
misfit. Finally, results are presented that show the interaction effects for various combinations of elastic 
edge restraint and initial misfit prestress. The values of the axial load P and buckling loads P cr are 
normalized by the linear-bifurcation buckling load of a geometrically perfect, quasi-isotropic panel with 
nominal dimensions and boundary conditions; that is, P cr ° = 15.78 kips (1 kip = 1000 lbs). Various 
displacement quantities presented herein are normalized by either the nominal panel thickness, t = 0.13 
in., the nominal panel length, L = 14.75 in., or the nominal panel width, W= 14.5 in. 

Elastic Circumferential Edge Restraint 

Results that show the effects of elastically restraining the circumferential movement of the loaded 
edges on the normalized buckling load are presented in figures 2-5. Specifically, load-shortening response 
curves for the geometrically perfect, quasi-isotropic panel with nominal dimensions and different degrees 
of circumferential restraint on the loaded edges are presented in figure 2. The corresponding curves for 
load versus normalized central out-of-plane displacement 8 It are shown in figure 3, where 8 = w(0, 0). 
Positive values of 8 It correspond to movement of the center of the panel outward, away from the center of 
curvature. Similarly, the corresponding curves for load versus normalized circumferential displacement 
vJW at one comer of the panel is shown in figure 4, where v a = v(L/2, - W/2 ). Negative values of vJW 
correspond to circumferential Poisson-like expansion that is normally associated with an axially 
compressed panel. Two curves are shown in each of figures 2-4 that correspond to values of the 
circumferential edge-restraint stiffness K v = 0 lb/in. (unrestrained) and K v = 10 7 lb/in., which is 
approximately 10.8 times the circumferential stiffness (A 22 ). A third curve is shown in each figure for the 
case of complete edge restraint, for which v = 0 at the loaded edges. The curves include the results for the 
prebuckling and postbuckling range of loading. The point at which buckling occurs for each case is 
marked on the curves with filled circles. Contour plots of the corresponding out-of-plane displacement 
field at the onset of buckling are shown in figures 5(a-c) for the three edge-restraint cases. The slight 
skewing that appears in the contour plots is associated with a relatively small degree of flexural 
anisotropy. 

The results in figures 2-4 indicate that the degree of circumferential edge restraint on the curved 
loaded edges of the compression-loaded panel has a significant effect on the response. For the case of 
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unrestrained edges (K v = 0), the load-shortening response curve shown in figure 2 is linear up to the 
buckling point at P/P cr ° = 0.98. Likewise, the corresponding center out-of-plane displacement shown in 
figure 3 is slightly inward prior to buckling, and the loaded edge of the panel undergoes circumferential 
expansion (figure 4), as expected for a circumferentially unrestrained edge. Near the buckling point, and 
in the postbuckling range of loading, the magnitude of the center out-of-plane displacement increases to 
many times the wall thickness of the panel. For the value of K v = 10 7 lb/in., the load-shortening behavior 
is different from the case of unrestrained edges. For this case, the load-shortening response curve (figure 
2) shows nonlinearity in the form of a slope change above P cr /P c ° = 1 , and it has a relatively high 
buckling load of P cr /P cr ° = 1.95. This nonlinearity in the response curve is caused by coupling of the 
circumferential and out-of-plane displacements near the loaded edges, which produce an outward bulge in 
the panel as the loading increases. The coupling of the displacement fields is indicated in figure 4 by the 
initial, slight circumferential expansion followed by circumferential contraction that is caused by local 
out-of-plane deformation. The outward bulge is evidenced by the corresponding monotonically increasing 
central out-of-plane displacement shown in figure 3 and the displacement contours shown in figure 5(b). 
The outward bulge acts to stiffen the panel, by increasing the moment of inertia in the central region, and 
raises the buckling resistance in the same manner as reported by Flui (ref. 4) for panels with a geometric 
imperfection in the shape of an outward bulge. For the case with fully restrained loaded edges (v = 0), the 
panel exhibits much less nonlinearity in its load-shortening response curve (figure 2), compared to the 
curve for K v = 10 7 lb/in. In addition, a less pronounced outward bulge (figures 3 and 5(c)) and a lower 
buckling load given by P cr /P c , ° = 1-52 is predicted for the panel with fully restrained loaded edges. 

Overall, the results in figures 2-5 show that prebuckling deformations retard the onset of buckling 
and result in increases in the buckling load of about 55% and 99%, with respect to the buckling load 
predicted for the panel with the loaded edges unrestrained, for the circumferential edge-restraint 
conditions defined by v = 0 and K v = 10 7 lb/in., respectively. Moreover, the results (figure 5) show that 
the prebuckling deformations are very sensitive to changes in the circumferential boundary condition. For 
example, the out-of-plane prebuckling displacements shown for the case of unrestrained edges (K v = 0) 
are characterized by a large, centrally located ellipse-shaped contours and an inward deformation pattern. 
For a value of K v = 10 7 lb/in., the predicted displacement contours are characterized by a hour-glass-like 
pattern, with maximum displacements on the order of 2.6 times the nominal panel wall thickness (figure 
5(b)). For v = 0, the displacement contours are characterized by a centrally located ellipse-shaped 
outward deformation pattern and the maximum displacement of about 1.5 times the nominal panel 
thickness (figure 5(c)). It is important to note that the results are contrary to intuition in that thay predict 
that the magnitude of the buckling loads and the severity of the nonlinear prebuckling deformations are 
not bounded by the corresponding results for the extreme boundary conditions defined by K v = 0 
(unrestrained) and v = 0 (completely restrained). 

The predicted effects of elastic circumferential edge restraint on the magnitude of the buckling 
load are presented in figure 6 as a function of the circumferential stiffness parameter K v . The buckling 
load values are marked in the figure by filled square symbols and can be placed into three distinct groups 
that were found to exhibit similar prebuckling behavioral characteristics. The first group of results is for 0 
< K v < 10 5 2 lb/in. Panels with edge-restraint stiffnesses in this group exhibit a linear prebuckling load- 
shortening response that is similar to that predicted for K v = 0, shown in figure 2. The normalized 
buckling load values for this group are all between 0.98 and 1.15. In addition, these panels exhibit 
prebuckling displacements similar to those shown in figures 3-5 for K v = 0. The second group of results 
are for panels with edge-restraint stiffnesses in the range 10 5 ' 2 lb/in. < K v < 10 7 2 lb/in. Panels in this group 
are predicted to exhibit some nonlinearity in their prebuckling response; that is, similar to the response 
characteristics shown in figures 2-5 for the panel with K v = 10 7 lb/in. Specifically, the prebuckling load- 
shortening response curves exhibit a nonlinear response above a load level of P/P cr ° = 0.95, and the 
buckling loads are on the order of two times the linear-bifurcation buckling load for the corresponding 
unrestrained panel. The third group of results is for panels with edge-restraint stiffnesses in the range K v > 
10 7 2 lb/in. The panels associated with this group are predicted to exhibit response characteristics similar 

6 



to those predicted for the fully restrained case (v = 0), as shown in figures 2-5. The portion of the curve 
shown in figure 6 for this group converges from above to the value of PJP C ° = 1-52, which is the value 
for the fully restrained panel. 

Elastic Rotational Edge Restraints 

Results that illustrate the effects of elastic rotational edge restraints on the buckling response of a 
compression-loaded curved panel are presented in figures 7-10. Specifically, load-shortening response 
curves for a compression-loaded panel with different degrees of rotational edge restraint on the straight, 
unloaded edges of the panel are presented in figure 7. The corresponding curves for load versus 
normalized central out-of-plane displacement 5 It are shown in figure 8, where 8 = w(0, 0). Three curves 
are shown in each of figures 7 and 8 that correspond to rotational edge restraint stiffnesses of K e = 0 
(equivalent to the classical simply supported boundary condition), K 0 = 1 0 3 in- lb, and the fully restrained 
case; that is, r x (x, - W/2 + 3/16 in.) = r x (x, W/2 - 3/16 in.) = 0. The buckling points are marked with the 
filled circles. 

The results in figures 7 and 8 indicate that the rotational edge restraint has a relatively small 
effect on the response and the buckling load of the panel. In particular, the load-shortening response 
curves are linear up to buckling and the buckling loads range from P/P cr ° = 0.98 to 1.08. The effective 
postbuckling stiffness increases slightly as the rotational constraint stiffness increases. This increase in 
postbuckling stiffness is caused by a reduction in the out-of-plane postbuckling deformations in the panel 
at a given load level, as shown in figure 8, and a subsequent increase in the effective load-carrying cross- 
section of the panel. The results in figure 8 indicate that the prebuckling center displacement of the panel 
is small and inward. However, near the buckling point, and in the postbuckling range of the loading, the 
center displacement increases to many times the wall thickness of the panel. For all cases, the predicted 
out-of-plane displacement contours at the onset of buckling are similar in character to that exhibited by a 
panel with nominal boundary conditions (K v = 0) and are characterized by a large ellipse-shaped inward 
deformation pattern (e.g., see figure 5(a)). 

A summary of the predicted effects of the rotational edge restraint on the buckling load of the 
panel is presented in figure 9. These results indicate that the buckling load increases very slightly from 
P cr /Per ~ 0.98 to 1.08, as K 0 increases from 0 to 10 9 in-lb, which is somewhat counterintuitive. This 
slight increase in the buckling load occurs in a transition zone for 10 2 in- lb < K e < 10 4 in- lb, in which the 
panel displacement response near the restrained edge transitions from the response for a simply supported 
edge to that for a clamped edge. More specifically, displacement traces for the panel with various 
rotational edge restraints are shown in figure 1 0 and illustrate the effects of the rotational restraint on the 
deformation response near the unloaded edge of the panel. These displacement traces are taken at the 
panel mid-length and span the panel arc from y = 0 to y = W/2. It is important to point out that the curves 
shown in this figure are not bounded by the curves that correspond to the classical simply supported and 
clamped boundary conditions because the curves correspond to different load levels incipient to buckling. 

For values of K e < 10 2 in- lb, the displacement response near the rotational restraint is similar to a 
free-rotation response that occurs with a classical simply supported boundary condition, as indicated in 
figure 10. For 10 2 in- lb < K e < 10 4 in- lb, the panel exhibits a significant reduction in the rotation near the 
unloaded edge and a reduction in the magnitude of the out-of-plane deformations at buckling, as 
compared to the simply supported boundary condition. For K 0 > 10 4 in- lb, the panel exhibits a similar 
displacement response to that of the fully restrained edge condition corresponding to r x = 0. 

Given the significant variation in the deformation response near the unloaded edge shown in 
figure 1 0, one might expect a much larger increase in the buckling load with an increase in the degree of 
rotational edge restraint. To address this expectation, and to understand better the panel behavior, a 
numerical parametric study of the effects of the panel (length-to-arc-width) aspect ratio, L/W, on the 
buckling response was conducted. Predicted buckling loads are plotted as a function of L/W for K 0 = 0 

7 



and r x = 0 in figure 1 1 . The curves shown in this figure are festoon curves, similar to the well-known 
corresponding festoon curves for buckling of flat, compression-loaded rectangular plates (e.g., see ref. 
11). The results indicate that the sensitivity of the panel response to variations in the degree of rotational 
restraint is clearly dependent on the panel aspect ratio. In addition, the results show that the particular 
panel configuration considered herein exhibits a relatively small difference between the buckling load for 
a panel with clamped unloaded edges and the buckling load of a panel with simply supported unloaded 
edges, at L/W = 1.02. Comparing these curves with the corresponding curves shown in ref. 11 for flat 
isotropic plates suggests that panel curvature amplifies the effect of clamping the loaded edges of the 
panels to the extent that the differences in buckling load associated with restraining the unloaded edges 
becomes less noticable. This behavior would account for the relatively small variation in buckling load 
with K 0 that is shown in figure 9. 

Initial Misfit Prestress 

Results that illustrate the effects of an initial misfit prestress (also referred to herein as an initial 
panel misfit) on the buckling response of the compression-loaded, quasi-isotropic curved panels 
considered herein are presented in figures 12-17. As discussed previously, the initial misfit prestress is 
caused by a mismatch between a 60-in-radius test fixture and a stress-free panel with a different initial 
radius R „ and is comprised of panel deformations and stresses that result from forcing the panel to 
conform to the 60-in-radius fixture. Displacement contours associated with the initial prestress of panels 
with initial radii R, equal to 52-in., 56-in., and 64-in. are presented in figure 12. Each panel configuration 
has two displacement contour plots associated with it. The first contour plot for each panel (figures 12(a)- 
1, 12(b)-l, and 12(c)-l) shows the predicted out-of-plane deformation that results from placement in the 
test fixture, as measured with respect to corresponding initial, undeformed geometry. The second contour 
plot for each panel (figures 12(a)-2, 12(b)-2, and 12(c)-2) shows the predicted deformed panel geometry, 
as measured with respect to a geometrically perfect 60-in-radius panel and is analogous to an initial 
geometric imperfection shape that is commonly used in shell buckling analyses. In general, the results 
indicate that the initial prestress deformations of the panels with R ; < 60 in. are outward bulges with 
magnitudes equal to 1-3% of the nominal wall thickness (see figures 12(a)-2, 12(b)-2). In contrast, panels 
with Ri > 60 in. exhibit initial prestress deformations that are inward bulges, as shown in figure 12(c)-2. 
In addition, initial inplane stress resultants and moment resultants were predicted for each case (not 
presented herein) and indicate that the stresses caused by forcing the panels into the test fixture are 
typically several orders of magnitude smaller than the corresponding stresses at the onset of buckling. 

Load-shortening response curves for panels with different amounts of prestress are presented in 
figure 13. The corresponding curves for load versus normalized central out-of-plane displacement 8/t are 
shown in figure 14, where 8 = w(0, 0). Similarly, the corresponding curves for load versus normalized 
circumferential displacement vJW at one comer of the panel are shown in figure 15, where v 0 = v{L/2, - 
W/2). The initial center displacements and circumferential displacements for P/P cr ° = 0 in figures 14 and 
15 are attributed to the initial prestress deformation response that occurs in the panels before the 
compression load is applied. Four curves are shown in each of figures 13-15 that correspond to panels 
with initial radii R, equal to 52 in., 56 in., 60 in. (no prestress), and 64 in., and illustrate the overall effects 
of an initial prestress on the compression response. The buckling points are marked with filled circles. 

The results in figures 13-15 indicate that the level of initial prestress in the panel has a significant 
effect on the prebuckling response. Panels with R, = 60 in. (no prestress) and 64 in. exhibit linear 
prebuckling load-shortening response curves and have normalized buckling loads of 0.98, and 0.90, 
respectively. Likewise, the corresponding center out-of-plane displacement shown in figure 14 is small 
and inward prior to buckling, with the displacements for the 60-in. panel being practically negligible, and 
the loaded edges undergo circumferential expansion, as shown in figure 15. The panel with R, = 56 in. 
exhibits a similar linear prebuckling load-shortening response and has normalized buckling loads of 1.17. 
However, the corresponding center out-of-plane displacement shown in figure 14 is small and outward 



prior to buckling. In addition, the load edge initially undergoes circumferential expansion followed by 
slight circumferential contraction at approximately P/P c ° = 1-02, as shown in figure 15. 

For a panel with an initial radius R f = 52 in., the behavior is significantly different. For this case, 
the panel exhibits a significant increase in the nonlinearity of the prebuckling response beyond a load 
level of P/P cr ° =1.2 and buckles at P/P c ° = 1.80. This nonlinearity in the response is caused by a 
coupling between the out-of-plane and circumferential displacements near the loaded edges of the panel, 
which produces an outward central displacement of the panel as loading increases. The strong coupling of 
the displacement fields is indicated in figure 15 by the initial, slight circumferential expansion followed 
by circumferential contraction that is caused by local bending. The outward prebuckling deformation is 
evidenced by the corresponding monotonically increasing central out-of-plane displacement shown in 
figure 14 and the displacement contours shown in figure 16(a). This response is similar in character to the 
response previously described for a geometrically perfect panel with a tangential, loaded-edge stiffness 
given by K v = 1.0 x 10 7 lb/in.(see Figs 2-5). Overall, the results for R, = 52 in. and 56 in. shown in figures 
13-16 indicate that the outward prebuckling deformations retard the onset of buckling and result in 
increases in the buckling load of about 84% with respect to the buckling load of the panel with no initial 
prestress. Moreover, the results in figure 16 indicate that the prebuckling displacement response can 
change greatly as the initial prestress in the panel changes. In particular, the results show that panels with 
Ri < 55 in. exhibit prebuckling displacements that are characterized by an hour-glass-like pattern with 
maximum displacements of approximately three times the nominal panel wall thickness, as shown in 
figure 16(a). This displacement pattern is very similar to that exhibited by the corresponding 
geometrically perfect panel with the circumferential edge restraint K v = 10 7 lb/in. that is shown in figure 
5(b). Panels with 55 in. < R , < 60 in. exhibit prebuckling deformations that are characterized by a large 
centrally located ellipse-shaped outward deformation pattern with maximum displacements of 
approximately 60% of the nominal panel wall thickness, as shown in figure 16(b). This displacement 
pattern is very similar to that of the corresponding geometrically perfect panel with the fully restrained 
circumferential edge condition, v = 0, shown in figure 5(c). Finally, panels with R, = 64 in. exhibit a 
prebuckling displacement response that is characterized by a large centrally located ellipse-shaped inward 
deformation pattern, as shown in figure 16(c), that is similar to that of a geometrically perfect panel with a 
circumferential edge restraint that ranges from 0 < K v < 10 3 lb/in., as shown in figure 5(a). 

The predicted effects of the initial, test- fixture induced prestress on the magnitude of the buckling 
load of the panel are summarized in figure 17. The buckling load values are marked in the figure by filled 
square symbols and can be divided into three groups; that is, R t < 55 in., 55 in. < R t < 60 in., and R t > 60 
in. In general, panels with R , > 60 in. exhibit a linear prebuckling load-shortening response and inward 
prebuckling displacements similar to those predicted for panels with R, = 60 in. and 64 in. that are shown 
in figures 13-16. The buckling-load values for this group range from P cr /P cr ° = 0.88 to 0.98. Panels with 
55 in. < R t < 60 in. exhibit linear prebuckling load-shortening response and outward prebuckling 
displacements similar to those predicted for panels with R, = 56 in., as shown in figures 13-16. The 
buckling load values for this group range from P cr /P c ° = 1.1 to 1.2. Finally, panels with R, < 55 in. 
exhibit a definite nonlinearity in their prebuckling response; that is, similar in character to the response 
shown in figures 13-16 for Rj = 52 inches. In particular, the prebuckling load-shortening response curves 
exhibit a large amount of nonlinearity above a load level of P/P er = 1.2, the prebuckling deformations are 
outward, and the buckling loads are on the order of two times the linear-bifurcation buckling load for the 
corresponding geometrically perfect panel with no initial prestress. 

Combined Effects of Initial Misfit Prestress and Elastic Edge Restraints 

Selected results illustrating the combined effects of an initial misfit prestress and elastic edge 
restraints on the compression response of a 60-in-radius, quasi-isotropic panel are presented in figures 18- 
21. More specifically, buckling loads for initial panel radii given by R t = 52 in., 56 in., 60 in., and 64 in. 
are presented in figures 18(a), 18(b), 18(c), and 18(d), respectively. Each figure presents normalized 
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buckling loads as a function of rotational restraint stiffness, K e , and contains seven curves that correspond 
to different degrees of circumferential restraint on the loaded edges of the panel. In general, the results 
indicate that there is a large amount of nonlinear coupling between the two elastic edge restraint 
conditions and the test-fixture-induced prestress. The results can be subdivided in to three major groups 
that correspond to various combinations of parameters K v , K e , and R, that are affecting the response. 
Group 1 is associated with panel responses that are affected by variations in K v only and correspond to 
panel configurations with K 0 > 10 4 in- lb and for all values of initial radius considered herein. The effect 
of K v on the panel response has already been presented herein and, as a result, will not be discussed in this 
section, since it was covered in one of the earlier sections. Group 2 is associated with panel responses 
that are affected by a combination of K v and initial prestress and corresponds to all panel configurations 
with K e < 10 2 in- lb. Group 3 is associated with panel responses that are affected by a combination of K v , 
K 0 , and initial prestress and include panel configurations with 10 2 in- lb < K e < 10 4 in- lb. The overall 
behavioral characteristics associated with the panels in Groups 2 and 3 are described next. 

Typical predicted effects of initial prestress and circumferential edge restraint on the buckling 
loads of panels in Group 2, panels that are affected by a combination of K v and initial prestress for K e < 
10 2 in- lb, are presented in figure 19. The results in this figure are for K e = 0. In general, the results 
indicate that the buckling loads of panels with R, > 56 in. exhibit a similar dependency on the 
circumferential edge restraint K v as geometrically perfect panels with the circumferential edge restraint 
only (see figures 2-5 and the related discussion). However, the prestress can affect the value of the 
circumferential edge restraint stiffness at which the panels transition from linear prebuckling behavior, 
associated with the relatively low buckling-load values of 0.98 to 1.15, to nonlinear prebuckling behavior, 
associated with the relatively high buckling-load values that range from 1.8 to 2.0. For panels with R, < 
54 in. and K v < 10 s lb/in., the response is primarily governed by the initial prestress and exhibits 
characteristics similar to those exhibited by the panels with only an initial prestress (see figures 12-17 and 
the related discussion). 

The response of the panels in Group 3 are affected by a combination of edge restraint conditions 
and initial prestress for values of 10 2 in- lb < K 0 < 10 4 in- lb. Within this group there are two subgroups; 
Group 3-1 for panels with K v < 10 5 lb/in. and Group 3-2 for panels with K v > 10 5 lb/in. The panels in 
Group 3-1 with R, = 52 in. exhibit nonlinear prebuckling behavior similar to that shown previously for a 
panel with K v = 10 7 lb/in. (see figures 2-5 and the corresponding discussion), and this behavior is 
characterized by large-magnitude, outward prebuckling deformations and buckling-load levels that 
approach 1.9 times the classical linear-bifurcation buckling load of the nominal panel, as shown in figure 
18(a). However, as R , increases, the panel response transitions to one that is characterized by a linear 
prebuckling response similar to that exhibited by a panel with nominal boundary conditions (see figures 
2-5 and the corresponding discussion for K v = 0), with relatively small inward prebuckling deformations 
and normalized buckling load levels between 0.90 and 1.08. Furthermore, the value of K 0 at which this 
transition from linear to nonlinear prebuckling behavior occurs is a function of the panel prestress and 
circumferential edge restraint K v , as indicated in figures 18(a) through 18(d). 

The buckling loads for panels in Group 3-2 (10 2 in-lb < K 0 < 10 4 in-lb and K v > 10 5 lb/in.) are 
influenced strongly by variations in the elastic edge restraints and to a lesser extent by the prestress, as 
shown figures 18(a) through 18(d). In particular, the curves for this class of panels exhibit a festoon-like 
behavior that is analogous to the curves associated with buckling of flat, compression-loaded, finite- 
length rectangular plates (for example, see Brush and Almroth, ref. 11). This festoon-like behavior is 
attributed to a distinct change in prebuckling and buckling response characteristics with an increase in the 
magnitude of the rotational edge restraint. For example, a panel with R, = 52 in., K v = 10 8 lb/in., and K 0 = 
10 2 in- lb exhibits response characteristics similar to those of a panel with v = 0 on the loaded boundaries 
(e.g., see figures 2-5). The response is characterized by a nonlinear prebuckling load-shortening curve, a 
normalized buckling-load value that ranges from P cr /P cr ° = 1.5 to 1.6, and a large centrally located, 
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ellipse-shaped deformation pattern at the onset of buckling (e.g., see figure 20(a)). As K e is increased to 
10 27 in- lb, the prebuckling response changes to become more like that of a panel with K v = 10 7 lb/in. 
(e.g., see figures 2-5). This response is characterized by a nonlinear prebuckling load-shortening curve, a 
normalized buckling-load value that ranges from P cr /P c ° = 1.80 to 1.98, and an hour-glass-shaped 
deformation pattern at the onset of buckling (e.g., see figure 20(b)). If K 0 is increased further to 10 4 in-lb, 
the prebuckling response changes once again. In particular, the buckling load is reduced from around 
1.98 to 1.5 and the magnitude of the maximum out-of-plane displacement decreases from 2.4 times to 1.0 
times the nominal panel wall thickness, as shown in figure 20(c). In addition, the character of the initial 
buckling displacements also changes. Specifically, for a panel with K 0 = 10 2 in- lb, buckling initiates near 
the loaded edges of the panel, as illustrated by the initial buckling displacement contours shown in figure 
21(a). 

These displacement contours illustrate, qualitatively, the change in the panel displacements at the 
onset of buckling and are determined by subtracting the initial buckling displacement field from the 
displacement field of previous solution step just prior to buckling and, as a result, are referred to herein as 
Aw contours. Displacement magnitudes are omitted from the plot because the displacement contours are 
intended to represent an infinitesimal change in the displacements in the neighborhood of the buckling 
point. For K 0 = 10 2 7 in- lb, the results indicate that buckling initiates near the unloaded boundary of the 
panel, as indicated by the initial buckling displacement contours shown in figure 21(b). In contrast, the 
initial-buckling displacements shown in figure 20(c) for the panel with K 0 = 10 4 in- lb show that buckling 
initiates in the center of the panel. 


Concluding Remarks 

A parametric study of the effects of test- fixture-induced initial prestress and elastic edge restraints 
on the prebuckling and buckling responses of a compression-loaded quasi-isotropic curved panel has been 
presented. Numerical results were obtained by using a geometrically nonlinear finite-element analysis 
code and high-fidelity models. Idealizations of the circumferential loaded-edge restraint and the 
rotational unloaded-edge restraint provided by the test fixture were included in the finite-element models 
to obtain insight into response-trend variations observed in previously conducted experiments. In 
addition, the effects of initial stress and deformations states, referred to herein as initial prestress, that 
result from a mismatch in radius of curvature between a test specimen and the test fixture were simulated 
analytically. 

Several of the effects of the test-fixture-induced initial prestress and elastic edge restraints have 
been presented, described, and discussed. In many cases, the results clearly indicate significant effects of 
the initial prestress and the two types of elastic edge restraint on the panel response. In particular, the 
results indicate that the initial prestress and elastic edge restraints can produce nonlinear prebuckling 
deformations in the panel that result in increases in the buckling resistance of the panel. In addition, the 
results indicate that various combinations of the initial prestress and elastic edge restraints can exhibit 
coupling and affect the overall nonlinear response of the panel in three different ways. Specifically, the 
predicted results can generally be categorized into three groups that exhibit distinct behavioral 
characteristics, depending on the relative degrees of the initial prestress and the elastic edge restraints. 
For some combinations of initial prestress and edge restraints, the panel exhibits a linear load-shortening 
response prior to buckling. The prebuckling deformations are negligible and the magnitude of the 
buckling load for these panels is typically very near the corresponding buckling load obtained from a 
linear-bifurcation buckling analysis. For other combinations of the initial prestress and edge restraints, the 
panel exhibits nonlinear prebuckling deformations and buckling loads that approach two times the 
corresponding linear bifurcation buckling load. For still other combinations of initial prestress and edge 
restraints, the panels exhibit smaller prebuckling nonlinearities and typically have buckling loads that are 
approximately 40% to 60% higher than the corresponding linear-bifurcation buckling load. 
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Overall, the results indicate, in a quantitative manner, how laminated-composite, curved-panel 
test specimens that are subjected to compression loads can exhibit significant sensitivity to the degree of 
circumferential and rotational edge restraint provided by a test fixture and the initial stress and 
deformation states that result from a small mismatch in the specimen and test-fixture geometries. 
Furthermore, since several combinations of initial prestress and edge restraints can result in similar 
response characteristics and buckling resistance, one must accurately characterize all pertinent aspects of 
the test fixture and panel specimen before each test in order to have the information needed to assess the 
influence of each factor on the panel response. 

In addition to the fundamental information presented herein that is relevant to structural testing, 
the results also provide insight into enhancements in buckling resistance that can be achieved by tailoring 
the support conditions for curved panels. The results should also be of particular interest to engineers and 
designers involved in simulating flight-hardware boundary conditions in structural verification and 
certification tests, and validating structural analysis tools. 
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Figure 2. Effects of circumferential loaded-edge restraint on the load-shortening response of a compression-loaded, 

60-in-radius, quasi-isotropic panel. 
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Figure 3. Effects of circumferential loaded-edge restraints on the out-of-plane center displacement of a 
compression-loaded, 60-in-radius, quasi-isotropic panel. 
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Figure 4. Effects of circumferential loaded-edge restraint on the circumferential of displacement of a compression- 

loaded, 60-in-radius, quasi-isotropic panel. 
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Figure 5. Effects of circumferential boundary stiffness Kv on the out-of-plane displacements at the onset of 
buckling of a 60-in-radius, quasi-isotropic panel (P cr ° = 15.78 kips). 



Figure 6. Effects of circumferential edge restraints on the buckling load of a compression-loaded, 60-in-radius, 

quasi-isotropic panel. 
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Figure 7. Effects of rotational unloaded-edge restraint on the load-shortening response of a compression-loaded, 60- 

in-radius, quasi-isotropic panel. 
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Figure 8. Effects of rotational unloaded-edge restraint on the out-of-plane center displacement of a compression- 

loaded, 60-in-radius, quasi-isotropic panel. 
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Figure 10. Effects of the rotational unloaded-edge restraint on the mid-length panel deformations of a compression- 

loaded, 60-in-radius, quasi-isotropic panel. 
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Figure 1 1 . Effects of panel aspect ratio L/W on the buckling load of panels with simply supported (K e = 0) and 
clamped (r x = 0) boundary conditions on the unloaded edges of a compression-loaded, 60-in-radius, quasi-isotropic 

panel. 
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(b-1) Rj= 56 in., referenced to the initial undeformed (b-2) Rj= 56 in., referenced to a geometrically perfect 
geometry 60-in-radius panel 
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Figure 12. 


Effects of initial radius Rj on the prestress displacements of a quasi-isotropic curved panel placed into a 

60-in-radius test fixture. 
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Figure 13. Effects of prestress on the load-shortening response of compression-loaded, quasi-isotropic curved 

panels in a 60-in-radius test fixture. 
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Figure 14. Effects of prestress on the out-of-plane center displacement of compression-loaded, quasi-isotropic 

curved panels in a 60-in-radius test fixture. 


25 



p/p 



Figure 15. Effects of prestress on the circumferential displacement response of compression-loaded, quasi-isotropic 

curved panels in a 60-in-radius test fixture. 
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Figure 16. Effects of prestress on the out-of-plane displacements of compression-loaded, quasi-isotropic curved 

panels in a 60-in-radius test fixture. 
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Figure 18. Effects of rotational boundary stiffness K 0 , circumferential boundary stiffness K v , and initial prestress on 
the buckling load of a compression-loaded, quasi-isotropic, curved panel. 
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Figure 19. Effects of initial prestress and circumferential stiffness K v on the buckling load of compression-loaded, 

quasi-isotropic, curved panels, K e = 0. 
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Figure. 20 Combined effects of circumferential boundary stiffness K v = 108 lb/in., rotational boundary stiffness Ke 
(in-lb), and prestress on the out-of-plane displacements at the onset of buckling of a 52-in. radius curved panel in a 

60-in. -radius test fixture. 



Figure 21. Combined effects of circumferential boundary stiffness K v = 108 lb/in., rotational boundary stiffness 
Ke (in-lb), and prestress the initial out-of-plane displacements at the onset of buckling of a 52-in-radius curved panel 

in a 60-in-radius test fixture. 
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